function [y1,y2] = dev_taylor(x, N)
    y1 = 1;
    y2 = 1;
    
    %Calcul des dev de Taylor y1 et y2
    for n = 1 : N
        y1 = y1 + x .^ n ./ factorial(n);
        y2 = y2 + (-1).^ n .*factorial(2.*n)./(1-2.*n) ./ (factorial(n)).^2./4.^n.*x.^n
    end
end


